Isotope effects and possible pairing mechanism in optimally doped cuprate 

superconductors 
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We have studied the oxygen-isotope effects on T c and in- plane penetration depth A a i>(0) in an 
optimally doped 3-layer cuprate Bii.6Pbo.4Sr2Ca2Cu30io+ s (T c ~ 107 K). We find a small oxygen- 
isotope effect on T c (ao = 0.019), and a substantial effect on A a b(0) (AA a f,(0)/A a (,(0) = 2.5±0.5%). 
The present results along with the previously observed isotope effects in single-layer and double- 
layer cuprates indicate that the isotope exponent ao in optimally doped cuprates is small while the 
isotope effect on the in-plane effective supercarrier mass is substantial and nearly independent of 
the number of the CuC>2 layers. A plausible pairing mechanism is proposed to explain the isotope 
effects, high-T c superconductivity and tunneling spectra in a consistent way. 



The pairing mechanism responsible for high-T c super- 
conductivity is still controversial. In conventional su- 
perconductors, a strong effect of changing ion mass M 
on the transition temperature T c implies that lattice vi- 
brations (phonons) play an important role in the micro- 
scopic mechanism of superconductivity. An isotope ex- 
ponent a (= —d\uT c /d\uM) of about 0.5 is consistent 
with the phonon-mediated BCS theory. A nearly zero 
oxygen-isotope effect (ao — 0.03) was earlier observed 
in a double-layer cuprate superconductor YBa2Cu307_ y 
which is optimally doped (T c ~ 92 K) Such a 

small isotope effect might suggest that phonons should 
not be important to the pairing mechanism. On the other 
hand, large oxygen-isotope shifts were later observed in 
several underdoped cuprate superconductors |§-§g]. Fur- 
ther, three indirect experiments have consistently demon- 
strated that the difference in the hole densities of the 16 O 
and 18 samples is smaller than 0.0002 per Cu site 
Moreover, a quantitative data analysis on the isotope- 
exchanged YBa2Cu306.94 suggested that there is a 
negligible oxygen-isotope effect on the supercarrier den- 
sity n s . 

Since muon-spin rotation experiments showed that 
T c is approximately proportional to n s jm** b in deeply 
underdoped cuprates (where rn** h is the in-plane effec- 
tive supercarrier mass), a large oxygen-isotope shift of 
T c observed in this doping regime should arise from a 
large oxygen- isotope effect on m** b . Indeed, several inde- 
pendent experiments |7[-p|,p^[ have consistently demon- 
strated that both the average supercarrier mass m** and 
the in-plane supercarrier mass m*£ strongly depend on 
the oxygen isotope mass in underdoped cuprates. Such 
an unconventional isotope effect suggests that there ex- 
ist polaronic charge carriers, which are condensed into 
supercarriers in the superconducting state. This appears 
to give a support to a theory of (bi)polaronic supercon- 
ductivity On the other hand, within this theory, 
it is difficult to explain a small isotope shift of T c and a 



large reduced energy gap (i.e., 2A(0)/fcsT' c > 6) observed 
in optimally doped cuprates where single-particle excita- 
tion gap vanishes above T c [Q. Therefore, an alternative 
theoretical approach is required to explain superconduc- 
tivity in optimally doped and overdoped cuprates. A 
possibly correct pairing mechanism should be able to con- 
sistently explain a small isotope shift of T c , a substantial 
isotope effect on the supercarrier mass ||,|l0| , and a large 
reduced energy gap. 

Here we report the observation of the oxygen-isotope 
effects on T c and in-plane penetration depth A a b(0) in 
a 3-layer cuprate Bii.gPbo.4Sr2Ca2Cu30io+ y (T c ~ 107 
K). We find a small oxygen-isotope effect on T c and a 
substantial effect on X a b(0). We propose a possible theo- 
retical model which is able to consistently explain these 
isotope effects and the large reduced energy gap. 

Samples of Bii.gPbo^S^C^CuaOio+y were prepared 
from high purity Bi 2 3 , PbO, SrC0 3 , CaC0 3 , and 
CuO. The samples were ground, pelletized and fired at 
865/855/845 °C for 37/40/44 h in air with two inter- 
mediate grindings. Two samples were prepared under 
nearly the same heat treatment. To ensure that the sam- 
ples have small grain size and enough porosity, they were 
reground thoroughly, pelletized and annealed in flowing 
oxygen at 600 °C for 10 h. 

Two pelletized samples were broken into halves (pro- 
ducing two sample pairs), and the halves were subject 
to the 16 and ls O isotope diffusion, which was con- 
ducted in two parallel quartz tubes separated by about 
I cm The diffusion was carried out for 68 h at 

600 °C and oxygen pressure of 0.8 bar for sample pair I, 
and for 40 h at 650 °C and oxygen pressure of 1 bar for 
sample pair II. The cooling rate was 30 °C/h. The oxy- 
gen isotope enrichments were determined from the weight 
changes of both 16 and 18 samples. The 18 samples 
had 85±5% 18 and 15±5% 16 0. 

The susceptibility was measured with a Quantum De- 
sign SQUID magnetometer. The field-cooled, measured- 
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on- warming susceptibility was measured in a field of 15 
Oe for sample pair I, and 10 Oe for sample pair II. The 
temperature measurements were performed with a plat- 
inum resistance thermometer (Lakeshore PT-111) placed 
in direct contact with the sample and driven by micropro- 
cessor controlled ac bridge in the SQUID magnetometer. 
The resolution is 2.5 mK and reproducibility is 10 mK at 
77 K after cycling to room temperature P,^p[. 
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FIG. 1. The susceptibility near T c for the 1B and ls O 
samples of Bii.6Pbo.4Sr2Ca2Cu30io+i/ (Bi-2223): (a) pair I; 
(b) pair II. 

In Fig. 1, we show the susceptibility near T c for the 
16 and 18 samples of Bii. 6 Pbo.4Sr 2 Ca2Cu30io+ y : (a) 
pair I; (b) pair II. In all the cases, the T c for the 18 
samples is 0.22±0.01 K lower than for the 16 samples. 
Extrapolating to 100% exchange, we calculate the iso- 
tope exponent a = -d\nT c J d\nM = 0.019±001. We 
also note that there is a well-dehned linear portion on the 
transition curve ~ 1 K below the diamagnetic onset tem- 
perature. It is evident that the slope of the linear portion 
(denoted by Pi) for the 18 samples is 5.0±1.0% smaller 
than for the 16 samples, that is, AP 1 /P 1 = -5.0±1.0%. 

For comparison, 
the results for single-layer Lai.ssSro.isCuC^ (LSCO) and 
double-layer YBa2Cu3C>6.94 (YBCO) compounds are re- 
produced in Fig. 2. It is clear that APi/Pi = -5.5±1.0% 
for LSCO and -6.8 ± 1.0% for YBCO. Comparing Fig. 2 
with Fig. 1, one can see that the isotope effect on Pi is 
nearly the same for all three compounds within the ex- 



perimental uncertainty. It is also interesting to note that 
the isotope effect on T c decreases monotonically with in- 
creasing T c . 

The observed oxygen-isotope effect on the slope Pi is 
caused by the dependence of the penetration depth on 
the oxygen mass. For nearly isotropic materials with 
A Q &(0) ~ A c (0) - A(0), it was shown that fuj: AP 1 /P 1 = 
-AT C /T C - 2AA(0)/A(0). For highly anisotropic mate- 
rials such as cuprates, a relation A C (T) >> A a b(T) > R 
(where R is the maximum particle size) holds near T c so 
that the diamagnetic signal is proportional to l/A^ fc (T) 
p5[ . Then one readily finds that 



APi/Px ~ -AT C /T C - 2AA afc (0)/A afc (0). 



(1) 



From Fig. 1 and Fig. 2, we obtain AA a b(0)/A a f,(0) = 
3.2±0.7% using Eq. 1. This is in remarkably good agree- 
ment with the recent muon spin rotation experiments 
on the oxygen-isotope exchanged YBa2Cu30g.96, which 
show that AA ab (0)/A ab (0) = 2.5±0.5% @. 
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FIG. 2. The susceptibility data near T c for the lb O 
and 18 samples of (a) Lai.ssSro.isCuCU (LSCO) and (b) 
YBa 2 Cu 3 06.94 (YBCO). After Refs. 

Since both n and n s are independent of the isotope 
mass as discussed above, the observed oxygen-isotope ef- 
fect on the in-plane penetration depth is caused by the 
isotope dependence of m* I . The substantial isotope effect 
on m*l may suggest that charge carriers in the optimally 
doped cuprates remain polaronic nature, and that those 
polaronic carriers are condensed into supercarriers in the 
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superconducting state. 

Now a question arises: why is the isotope effect on m*£ 
substantial while the isotope shift of T c is very small in 
these optimally doped cuprates? In order to answer this 
question, we should first find out which phonon modes 
are strongly coupled to doped holes in these materi- 
als. Inelastic neutron scattering experiments jl7],[l8| show 
that the Cu-0 bond stretching mode in the Cu0 2 planes 
are strongly coupled to the doped holes. The average fre- 
quency of this mode is about 75 meV in Lai.gsSro.isCuC^ 
0, and about 60 meV in YBa 2 Cu 3 6 .92 @. Such 
a strong electron-phonon coupling should be also man- 
ifested in tunneling spectra, as it is the case in con- 
ventional superconductors. Although it is difficult to 
obtain reliable tunneling spectra for cuprates due to a 
short coherent length, there are two high-quality tun- 
neling spectra for slightly overdoped YBCO jl9| and 
Bi 2 Sr 2 CaCu 2 8+a (BSCCO) @. 
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FIG. 3. The normalized conductance data for a scanning 
tunneling microscopy (STM) tunnel junction on a slightly 
overdoped YBa2Cu307_ y (YBCO) crystal with a Pt-Ir tip 
at 4.2 K. After [[TJ). 

In Fig. 3, we show normalized conductance data for a 
scanning tunneling microscopy (STM) on a slightly over- 
doped YBCO crystal with a Pt-Ir tip at 4.2 K Jig . The 
crystal has T c ~ 90 K with ~ 1 K transition width jl9) . 
Since the presence of oxygen vacancies in the CuO chains 
can lead to residual density of states and to zero-bias 
conductance in the superconducting state, a negligible 
zero-bias conductance in the spectrum suggests that the 
spectrum represents the intrinsic density of states con- 
tributed only from the Cu0 2 planes. It is striking that 
the strong coupling features similar to that in conven- 
tional superconductors can be clearly seen in the spec- 
trum (as indicated by the arrows). The strong-coupling 
features correspond to a strong coupling between charge 
carriers and the phonon mode with u> p h — 19 meV. The 
phonon density of states in an optimally doped YBCO 



also reveals a large peak at about 20 meV |H) . 

In Fig. 4, we plot the electron-phonon spectral density 
a 2 F(u>) for an optimally-doped BSCCO crystal, which 
was extracted from an SIS break-junction spectrum [po[ . 
A strong coupling feature at an energy of about 20 meV 
is clearly seen. This feature also corresponds to the large 
peak in the phonon density of states at about 20 meV 
(see open circles). In addition to a strong coupling fea- 
ture at 20 meV, there is another strong coupling feature 
at about 73 meV, which corresponds to the phonon en- 
ergy of the Cu-0 bond stretching mode discussed above. 
Therefore, these tunneling spectra (Fig. 3 and 4) con- 
sistently suggest that both the low-energy phonon mode 
at about 20 meV and the high-energy phonon mode at 
about 73 meV are strongly coupled to conduction elec- 
trons in these double-layer compounds. 
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FIG. 4. The electron-phonon spectral density a 2 F(uj) for 
an optimally doped Bi2Sr2CaCu2 08+;, (BSCCO) crystal, 
which was deduced from an SIS break-junction spectrum [^o) . 

A theoretical approach to a strong electron-phonon 
coupling system depends not only on the adiabatic ra- 
tio uj/Ep but also on the coupling strengths with dif- 
ferent phonon modes. When the Fermi energy Ep is 
smaller than these strong-coupling phonon energies, the 
phonon-induced effective interaction between carriers is 
nonretarded so that the real-space pairing (g-R^, inter- 
site bipolaron formation) becomes possible [|13|,|22|. This 
should be the case for the doping level x < 0.10 in 
Lai_ 2; Sr 2 ;Cu04 In contrast, Ep in doped oxygen- 
hole bands may lie in between 20 and 73 meV in the 
optimally doped and overdoped regimes, so that the 
pairing interaction becomes retarded for the low-energy 
phonons, and remains nonretarded for the high-energy 
phonons. The retarded electron-phonon coupling for the 
low-energy phonons could be treated within the Migdal 
approximation, while the nonretarded electron-phonon 
coupling for the high-energy phonons should be mod- 
eled separately within the polaron theory. This theoret- 
ical approach has been successfully applied to fullerenes 
p3|. The strong coupling between doped holes and the 
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high-energy phonons leads to a polaronic mass enhance- 
ment and to an attractive nonretarded potential between 
doped holes. Effectively, the polaronic holes could then 
form k-space Cooper pairs by interacting with the low- 
energy phonons. The problem could thus be solved 
within Eliashberg equations with an effective electron- 
phonon spectral density for the low-energy phonons and 
a negative Coulomb pseudopotential produced by the 
high-energy phonons and other high-energy bosonic ex- 
citations of purely electronic origin (e.g., spin fluctua- 
tions, excitons, and plasmons). Within this simplified ap- 
proach, the effective electron-phonon coupling constant 
A ep for the low-energy phonons is enhanced by a factor of 
f p — exp(g 2 ). Here g 2 = A/ujh, A is a constant, and u>h 
is the frequency of the high-energy phonon mode. The ef- 
fective Coulomb pseudopotential is negative and also 
proportional to f p . 

9 

from the mid-infrared optical conductivity which exhibits 
a maximum at E m ~ 0.12 eV (26). With E m = 0.12 
eV, huj H = 75 meV, we find g 2 = E m /(2Kuj H ) = 0.8, 
leading to f p = 2.2. From the spectral density shown 
in Fig. 4, we can extract the effective electron-phonon 
coupling constant X ep for the low-energy phonon mode, 
that is, A ep — 2.6. If there were no polaronic mass en- 
hancement due to the high-energy phonons, the coupling 
constant contributed from the low-energy phonons would 
be about 1.2. With fi* = 0.1 and X ep — 1.2, we calculate 
T c = 18 K according to a T c formula 



For slightly overdoped BSCCO, g 2 can be evaluated 



k B T c = 0.25W< uj 2 >[exp(2/A e// ) - 1] 



-1/2 



(2) 



where 



(3) 



^eff = (A ep - AO/[l + 2 A** + A ep /X*t(A e p)], 

The function t(X ep ) is plotted in Fig. 2 of Ref. ^7|. In the 
present case, ftV < lj 2 > is contributed only from the low- 
energy phonons and equal to 20 meV. Therefore, without 
the high-energy phonons, T c would not be higher than 20 
K. The high-energy phonons not only enhance A ep by a 
factor of 2.2, but also reduce fi* substantially [^3). It has 
recently been shown that fi* in cuprates becomes nega- 
tive (i.e., [i* ~ —0.05) due to the presence of low-energy 
electronic collective modes (acoustic plasmons) in layered 
conductors j23]. Since the high-energy phonon mode re- 
duces \i* further, it is likely that the value of fi* should 
be in the range of —0.1 to —0.2. If we take fi* = —0.15, 
A e p = 2.6 (see above), we can get T c — 95 K. This leads to 
k B T c /(hy/< uj 2 >) = 0.41, and 2A(0)/k B T c g 7 ac cord- 
ing to the known relation between ksTc/ '(hy/ < ui 2 >) and 
2A(0) /ksT c Ea] . The calculated reduced energy gap is in 
good agreement with experiment. We would like to men- 
tion that the present calculation is valid for an isotropic 
s-wave gap, which should be the case for the polaronic 
oxygen-hole bands near (0,±7r) and (±7r, 0) regions. In 



addition to the polaronic oxygen holes, there are electron- 
like free carriers which could condense into supercarriers 
through interband scattering. The pairing symmetry of 
the two-carrier system may be an extended s-wave with 
eight line nodes |2S| , 

Now we can calculate the total isotope exponent a us- 
ing Eqs. H and ||, and the relations: A ep oc f p , pb* cx / p , 
t(X ep ) ~ 1.8/Aep. The calculated total isotope exponent 
is a ~ 0. The nearly zero isotope exponent is due to the 
fact that the isotope dependencies of A ep and fi* (arising 
from the polaronic effect) cancel out the isotope effect on 
the prefactor of Eq. ||. Moreover, the formation of pola- 
ronic Cooper pairs naturally explains the sizable isotope 
effect on the supercarrier mass. 

In summary, we report the observation of the oxygen- 
isotope 

effects in the 3-layer cuprate Bii.6Pbo.4Sr2Ca2Cu30io+ tf . 
The present results along with the previously observed 
isotope effects in single-layer and double-layer cuprate 
superconductors indicate that the isotope effect on T c 
in optimally doped cuprates is small while the isotope 
effect on the effective supercarrier mass is substantial. 
These isotope effects and tunneling spectra observed in 
optimally doped cuprates can be consistently explained 
within a scenario where polaronic oxygen holes are bound 
into Cooper pairs. 
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